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Resolution of both approaching and receding ejecta in the galactic microquasars makes possible 
to measure the flux ratio S'a/ S-c of twin ejecta, which contains an important information about the 
pi ^. nature of the jets. We show that the flux ratio S^jS-c = 8 ± 1 observed from GRS 1915+105 during 

' the prominent March/ April 1994 radio flare can be explained in terms of relativistic motion of 

discrete radio clouds, if one assumes that the twin ejecta are similar, but not completely identical, 
i.e. allowing for some asymmetry between the plasmoids in their speeds of propagation and/or 
luminosities. The recoil momentum due to asymmetrical ejection of the pair of plasmoids could be 
comparable with the momentum accumulated in the inner accretion disk. We suggest a possible 
explanation of the observed anticorrelation between the X-ray and radio flares as the result of 
Q> , drastic structural changes in the inner disk caused by production of powerful jets. The delay 

I between the times of decline of X-rays and appearance of strong radio flares is explained by time 

' needed for expansion of radio clouds to become optically thin. 

■ Subject headings: stars: individual (GRS 1915-1-105) - stars: flare - galaxies: jets - radio contin- 

— I uum: general - accretion disks 

O ■ 1 Introduction 

. Recent discovery of the galactic superluminal jet sources GRS 1915-1-105 (Mirabel & Rodriguez 1994, 

hereafter MR94) and GRO J1655-40 (Tingay et al. 1995; Hjellming & Rupen 1995) provides unique 
^ , possibility for deep study of the phenomenon of relativistic jets observed elsevi^here in the Universe. 

These sources, called microquasars, represent, to a large extent, scaled down analogs of AGNs (Mirabel 
&: Rodriguez 1995). Being, however, much closer to us than AGNs, the microquasars enable radio 
monitoring of both approaching and receding relativistic ejecta, and importantly, in short timescales. 

GRS 1915-1-105 is the first source of apparent superluminal jets where the receding jet has been 
also detected, and time evolution of radio fluxes of both components of the prominent 19 March 1994 
flare has been traced on timescales of days beyond 30 April (MR94). In particular, this has allowed 
Mirabel & Rodriguez to determine the real speed f3 ^ 0.92 and the angle of propagation 9 « 70° of 
the ejecta. At an estimated distance D = 12.5 it 1.5 kpc, the observed proper motion of radio sources 
correspond to apparent velocities (in units of speed of light) fa = 1.25 ± 0.15 and Vr = 0.65 ± 0.08 for 
the approaching and receding ejecta, respectively (MR94). 

Detection of pairs of jets makes possible to put forward an important question concerning the origin 
of the radio images observed. Namely, interpretation of the ratio of radio fluxes detected from the 
twin jets in GRS 1915-1-105 actually addresses an important issue, whether the moving radio images 
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are the radio 'echos' of relativistic shocks propagating in the jets, or they correspond to real motion 
of discrete plasmoids. 

The ratio of flux densities, S{i') oc measured at equal angular separations of the identical pair 
of jets from the core is 

Sa _ ( I + /3cOsg y+" 

S, \l - f3 cose) ' ^ ^ 

The index k is equal to = 3, if the fluxes are produced in the moving discrete radio clouds, but 
k = 2 for the brightness ratio of continuous stationary jets (e.g. see Lind &: Blandford 1985). The 
difference in the index k is due to relativistic volume contraction effect (see Bodo &: Ghiscllini 1995). 
For P = 0.92, e = 70°, and measured radio spectral index a = 0.84 ± 0.03 (MR94), Eq.(l) results 
in the flux ratio Sa,/Sr 2± 12 to be expected from discrete sources, whereas the measured flux ratio 
'S'a/'S'r = 8 =b 1 (MR94) is closer to -SaZ-Sr ~ 6 which is formally expected in the case of stationary jets. 
Meanwhile, the apparent motion of discrete radio condensations does not support the interpretation 
of the observed flux ratio in terms of emission of stationary jets. A possible explanation of this 
discrepancy between the observed and expected (for discrete plasmoids) flux ratios of the jets in GRS 
1915+105 was given by Bodo &: Ghisellini (1995), who suggested that the real speed of radioemitting 
plasma (fluid) in the jets is different from the speed /3 0.92 of radio patterns attributed to propagation 
of shocks in the fluid. 

In this paper we propose another interpretation of the measured flux ratio, in terms of speeds of 
the ejecta coinciding with the speeds /3 ~ 0.9 of radio patterns, assuming that the twin ejecta are 
similar, but not absolutely identical. 



2 Pair of asymmetrical jets 

When both components of relativistic jets can be detected, then in addition to the equation for the 
angular speed //a of the approaching jet, two more equations become available, namely, the one for 
the angular speed /Xr of the receding component, as well as the equation for the flux ratio of the two 
components. Allowing for the parameters of the counter ejecta to be different, the first two equations 
read: 

_ /3a sin 6>a _c 

1 - /3a COS 0a D 

_ fi- sin Or c_ 

l + (5rCose, D ' ^"^^ 
where /3a and /3r are the speeds, and 0^ are the angles of propagation of approaching and receding 
ejecta. Note that for convenience in Eq.(3) we have substituted 9^ 180° — Ox- 

The equation for the flux ratio Sa,/ Sr for an asymmetrical pair of ejecta is found, using the relation 
S{i') = 6^^'^ S'(y) (e.g., Lind & Blandford 1985) between the apparent and intrinsic energy fluxes of 
a radio cloud moving with Dopplcr factor 5 = \/\ — /3^/(l — (3 cos 9). Then the flux ratio at equal 
intrinsic times is: 

5a ^ Tr \ 3+" / 1 + (3, cos 9, \ L', 



S, \T^) \l-l3^cos9j 
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where Fa^r = (1 ~ /^a.r)"^^^ the Lorcntz factors of the bulk motion, and L^r are the intrinsic 
luminosities of the clouds. An asymmetry between the ejecta would imply at least one of the following 
inequalities: (a) the luminosities of the approaching and receding ejecta in their rest frame are 
different, i.e. / ; (b) the ejecta propagate not strictly in the opposite directions, i.e. 9a, Oi] 
(c) the speeds of the ejecta are different, Pa. Pr ■ In principle, one may assume also that the spectra 
of the ejecta are described by different power-law indexes a. However, since the same a ^ 0.84 was 
observed from both radio clouds of the 19 March 1994 flare, we do not consider here such a possibility. 

Eq.(4) corresponds to the flux ratio of two radio clouds at equal intrinsic times, t'.^ = t[. = t'. 
Meanwhile, the directly measurable quantity is the flux ratio at equal angular separations (/>a = 0r = 0- 
The relation between these two quantities, (Sa/Sr)t' and {Sa/Sr)cf,, can be found taking into account 
that in the observer's frame the ratio of the times corresponding to equal angular separation of the 
sources from the core is (ta/^r)^ = IJ-i/fJ-aj while the ratio of apparent times ta,r = t' /5a,r corresponding 
to the same intrinsic time t' is {ta/tr)t' = Sr/Sa- Using Eqs. (2) and (3), we find that (ta/ir)</) = 
Hta/tr)t' where 

_ Tr /3r sin 9r , . 

~ ra/3aSin^a ' ^ ' 

Approximating the time evolution (decline) of the flare as 6'a,r oc t~P, as it was for the 19 March 1994 
flare of GRS 1915+105 (with p ~ 1.3, Mirabel Sz Rodriguez 1995), the relation between the flux ratios 
at equal intrinsic times and equal angular separations is found: 

(I). ^ (I). ■ 

For GRS 1915+105 calculations below result in A ~ 0.85, therefore the measured {Sa/Sr),j> ~ 8 
corresponds to {Sa/ Sr)t' ~ 6.5. Further on we will omit subscript (i'), implying under {Sa/Sr) the 
flux ratio at equal intrinsic times, if not mentioned otherwise. 

Returning now to Eqs. (2)-(4), note that for identical jets, when f3a = Pr = /3, ^a = = 9i Eq.(5) 
results in A = 1, therefore in this particular case the equal angular separations correspond to equal 
intrinsic times, as expected, so Sa/ = (Sa/Sj:)^. Taking into account also that L'^ = L'^. , Eq.(4) 
comes to Eq.(l) with k = 3. Then Eqs. (l)-(3) make formally a system of 3 equations for 3 variables: 
/?, 0, D. However, these equations are not independent when the power-law index k is fixed, since they 
predict the ratio -Sa/^r = (/Xa//Xr)'^"'"". Note that due to high value of the index A; + a ~ 4, the predicted 
flux ratio is very sensitive to variations of the angular velocities /ia,r of the twin ejecta. In particular, 
in the case of /c = 3 and a = 0.84, small error limits of the observed values of = 17 ± 0.4mas/d 
(milliarcsec/day) and /Ur = 9.0 ± O.lmas/d (MR94) result in a rather broad range of the flux ratios 
predicted for identical plasmoids, {Sa/S^)^ = Sa/S^ = 13.1^\'^ , which is, nevertheless, beyond the 
observed (^a/Sr)^ = 8 + 1. 

However, this discrepancy disappears when an asymmetry between the ejecta is allowed. Consider 
separately each of the 3 basic types of asymmetries mentioned above. 

A. f3a = Pt = P and 6^ = 0^ = 6 , but La 7^ Lj,. In this case the solution to Eqs. (2) and (3) for GRS 
1915+105 is given in MR94: (3cos9 = 0.323, and for the distance D = 12.5 kpc the true speeds of the 
ejecta /3 0.92 and 9 70°. Then from Eq.(4) follows that in order to explain the ratio (>S'a/S'r)</) ~ 8, 
the ratio of intrinsic luminosities should be L'^/L'^ ~ 1.6. 
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B. /3a = /5r = /? and L'^ = L[. , but Og^ ^ 6^. Introducing dimensionless parameters Vg^ = fi^D/c and 
Vj: = fij:D/c, which correspond to the apparent speeds (in units of speed of Ught) of the ejecta for a 
given distance D, the solution to Eqs. (2)-(4) can be written as: 

2 (s — 1) Va_ 

sin Or = b sin , (8) 

Va 

P = TT > (9) 

sm + Va. COS 9^ 

where s = (S'a/S'jVCs+a) ^ and 6 = svjva, = snj/ia,. Eqs. (7)-(9) define 6'a, 0^ (in the range < 180°) 
and P for arbitrary /Xa, /^r and s (note that if 6 > 1, one should care for the condition sin^?;. < 1). 
However, starting from some distances these equations result in /3 > 1, which becomes formally possible 
after reduction of Fa and Tj. in Eq.(4). For the values of /^a smd fij. observed for the March 19 event, 
and Sa/Sj- < 7, the solutions with /3 < 1 require D < 10.7 kpc. For example, assuming D = lOkpc 
and Sa/Sr = 7, Eqs. (7)-(9) result in 6'a = 83.7°, 9^ = 57.7° and p = 0.92. From Eq.(8) follows that 
symmetrical solutions, with 9^, = O^, exist only if 6 = 1, which is not the case in GRS 1915+105. 

C. La = and 9^ = 9r = 9 , but Pa^ Pr- In this case, due to very strong dependence of the flux 
ratio Sa/Sr on the ratio Fa/Fr (see Eq.(4)), the discrepancy between the observed and 'expected' flux 
ratios by factor of (1.5 — 2) could be easily explained with only '-^ (10 — 20) % difference in Lorentz 
factors of the oppositely moving plasmoids. Dividing Eq.(2) to Eq.(3), and substituting the resulting 
equation into Eq.(4), we find the relation between /3r and Pa, and afterwards Eqs. (2) and (3) can be 
resolved with respect to Pa and x = Pa cos 9: 

P. = ^|x^ + vl{l-x? , (11) 

Parameters s and b are the same as in Eqs.(7)-(9). The angle 9 is readily found after calculations of 
X and /Ja- 
in Fig.l we show the dependence of the speeds of propagation, (3^ and Pr, of the approaching and 
receding plasmoids on the distance D, calculated from Eqs. (10)-(12) for //a = 17mas/d, = 9mas/d, 
and fixed flux ratio at equal intrinsic times Sa/S^ = 7. For the distance D = 12.5 kpc the speeds 
Pa = 0.926 and /?r = 0.902, and the angle 9 = 70.2° are found. Prom Eq.(6) follows that in this case 
the flux ratio S'a/5'r = 7 corresponds to {Sa/Sr)cf, — 8.6. Remarkably, the speeds of both ejecta are 
close to the value of /3 = 0.92 found in MR94 assuming j3a = Pt- In Pig.2 we show the speeds Pa and 
Pr and relevant Lorentz factors Fa and Fr of counter ejecta, calculated for the same /ia,r and a as in 
Fig.l, but for different fiux ratios in a broad range 5 < Sa/ < 15. Note that small power-law index 
1/(3 + a) = 0.26 in the parameter s = (•S'a/S'r)^^^^"''"^ significantly reduces the impact of uncertainties 
in the flux ratio Sa/S^ on calculated speeds of propagations of the ejecta. 
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3 Discussion 



The key assumption of this paper consists in allowance for an asymmetry between the pair of plasmoids, 
which implies at least one of the enequalities: 7^ L^, Og, 7^ ^r, or /?a 7^ /?r- Assuming that only the 
angles of propagation are different, Og, ^ 9^, one can explain the observed flux ratio {Sa/Sr)(f, = 8 ± 1 
if GRS 1915+105 would be at a distance D ~ lOkpc, and requiring essential difference (> 25°) 
between and 9j-. Note, however, that in this case the distance to GRS 1915+105 does not agree 
with estimated kinematic distance D = 12.5 ± 1.5 kpc (MR94). Besides, the radio observations of 
the March 19 outburst clearly show the motion of two sources in strictly opposite directions on the 
sky, within uncertainties of only few degrees (MR94) . Therefore the difference in the angles 9^ and 9r 
could be only in the plane of view. 

More probable seem to be realization of the other two types of asymmetries. If the asymmetry 
is attributed mainly to different intrinsic (i.e. in the rest frames of the clouds) luminosities L'^ ^ L'^., 
then the receding cloud should be significantly (> 50 %) more luminous than the approaching one. In 
principle, such a difference in the luminosities of approaching and receding clouds could be attributed 
to different (by ~ 30%) magnetic fields B, and/or to different amount of relativistic electrons in the 
ejecta. At last, the observed flux ratio can be explained with only a small difference in the speeds 
of propagation, Pg, ^ P,-, corresponding to ~ 10% difference between the Lorentz factors Fa and Fr 
of the bulk motion of the plasmoids. Due to very strong dependence of Sg,/Sr on the ratio Fa/Fr, a 
small asymmetry in the speeds of propagation of the jets can be considered as the main reason for 
the discrepancy between the measured and 'expected' flux ratios in GRS 1915+105. Note, however, 
that there is no particular reason to believe that only the velocities of the ejecta are different, and it 
is quite possible that the intrinsic luminosities of the clouds can be also somewhat different, L'^ ~ L'^ 
(but 9g, ^ 9r). In that case Eqs. (10)-(12) can be used for calculations of the jet parameters, if one 
substitutes parameter s = (S'a/S'r)V(3+«) there with = {Se^L'JSr L'J^^(^+''\ 

The asymmetry in the luminosities and/or speeds of propagation of the ejecta could be connected 
both with somewhat different conditions in the external medium at distances < 10^^ — 10^^ cm from 
the core (reached by the ejecta in a few days) or with an asymmetry of the bidirectional ejection itself. 
In the absence of a firmly established theoretical model for development of the two-sided jets, there 
is no particular reason to eliminate any of these two options. For example, different densities of the 
external medium might result in different speeds of propagation, and/or different speeds of expansion 
of radio clouds leading to different magnetic fields in the clouds at equal intrinsic times. An asymmetry 
of the bidirectional ejection (which is obvious for GRO J1655-44, see Hjellming & Rupen 1995), e.g. 
~ 10% difference in the masses or acceleration rates of the ejecta, could also result in different speeds 
of propagation and/or luminosities of the plasmoids. 

Remarkably, intrinsic asymmetry of bidirectional ejection process could imply an interesting inter- 
pretation of the anticorrelation between the X-ray and radio flares observed from both microquasars 
GRS 1925+105 (Foster et al. 1996; Harmon et al. 1997) and GRO J1655-44 (Tingay et al. 1995; 
Harmon et al. 1995; Hjellming & Rupen 1995), which consists in appearance of strong radio outbursts 
generally after essential decline of the preceding X-ray flares, usually with some delay up to few days 
between these events. Tingay et al. (1995) and Meier (1996) interpreted this effect as an indication of 
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suppression of jet production mechanisms in the regime of supercritical accretion responsible for strong 
X-ray flares, while ejection would occur only after slowing down of the accretion rate to subcritical 
regime. Possible asymmetry of the jets allows us to suggest an alternative (in fact, to a large extent 
opposite) scenario where the jet production itself may result in temporary drop of the accretion from 
supercritical to subcritical regimes. 

Indeed, internal energy in the radio clouds, and correspondingly, minimum kinetic energy of the 
ejection for strong outbursts in GRS 1915+105 is estimated as ~ lO^^erg (Liang &; Li 1995; 
Meier 1996), and even more (MR94). Then the momentum transferred to relativistic ejecta can be 
estimated as Pa,r ~ W/c > 3 x 10^^ gem /s. Asymmetrical ejection generally would require, just from 
momentum conservation law, that significant recoil momentum AP = |Pa — Pr| is transferred to the 
core of ejection, which is the third object in the interacting system "two jets + core". Since most 
probable site responsible for relativistic ejection is the inner part of the accretion disk (e.g., Blandford 
& Payne 1982; Begelman et al. 1984), it is interesting to compare AP with the integrated value of 
the specific (per unit volume) momentum pv = p\v\ of the gas orbiting in the inner accretion disk 
at radfi R' < R: Pd{R) = J pvd^R' (here p is the mass density in the accretion flow). Using the 
Shakura & Sunyaev (1973) accretion disk model, this integrated 'absolute momentum' is estimated 
(in their a-region) as Pa — 4 x 10^^ a^g m^^ gem /s, where ass < 1 is the viscosity parameter, 
and r = R/R^ is the radius in the units of i?g = 3 (M/Mq) km. Then, even for ass ~ 0.01, the black 
hole mass m = M/Mq ~ 10, and accretion rate about the critical, m ~ 1, a small recoil momentum 
AP ~ 0.1Pa,r being applied to the jet production region (approximately, perpendicular to the disk 
plane) would be enough to cause significant changes, or even destruction, of the inner disk up to radii 
r* ^ 100 estimated from condition AP ~ Pd(r*), and resulting in a temporary reduction/termination 
of the fuel supply into this region responsible for thermal X-rays. 

Thus, in the framework of this scenario, the onset of subcritical/supercritical accretion would 
correspond to the active state, with X-ray flares which could proceed without significant ejection 
events and observable radio flares (see Foster et al. 1996; Tavani et al. 1996). Powerful ejection of 
radio emitting material may be accompanied by significant destruction of the inner accretion zone, 
leading to a strong decline of the X-ray fluxes simultaneously with production of relativistic ejecta. 
The ejection event, however, will be not accompanied by the simultaneous increase of the radio fluxes, 
since for an appreciable time the ejecta remain too compact to be transparent at the radio frequencies. 
Thus, one has to expect a delay, up to a few days, between the instant of ejection and the appearance 
of strong radio flares, which time is needed for expansion of the plasmoids, with Ve^cp ~ 0.1c, to the 
radii Pel > 5 x 10^^ cm, as follows from the optical transparency of the radio plasmoids with respect 
to synchrotron self-absorption (Atoyan &; Aharonian 1997). Note that the radio clouds, moving with 
the speeds ~ c, become transparent at distances P ^ 10 Pd > 5 x 10^^ cm, i.e. far away from the 
binary system. Depending on the timescale needed for the accretion disk to recover after ejection 
event, the maximum of radio flares could appear just in the dips between X-ray flares (see Harmon 
et al. 1997). Rough estimate of the duration At of such dips as the dynamical time of accretion from 

7/2 

radii r^, results in At oc rj , which could vary on time scales from minutes to days, depending on the 
spatial scales of the inner disk affected by ejection. 
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Figure 1: The speeds of propagation of the approaching (soUd Hne) and receding (dashed hnc) 
asymmetrical pair of plasmoids calculated from Eqs.(10)-(12) assuming different distances to GRS 
1915+105. 
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Figure 2: The speeds (bottom panel) and the Lorentz faetors (top panel) of the bulk motion of the 
pair of plasmoids calculated for a fixed distance D = 12.5 kpc and angular velocities //a,r observed 
from GRS 1915+105 (MR94), but assuming different flux ratios Sa,/Sr. 
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